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We analyze some signatures of neutral Higgs bosons, produced in high energy pro- 
ton (anti)proton collisions, which decay primarily to top quarks. This channel is 
^ ■ clouded by a large irreducible it background. We investigate the it invariant-mass 

distribution in the lepton+jets decay mode of it and a it spin-spin correlation in 
the dilepton decay mode. At the LHC such Higgs bosons can be detected in the it 
invariant-mass spectrum, and the spin-spin correlation studied by us is also sensitive 
to a Higgs boson if this particle has a sizeable pseudoscalar component. 
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1 Introduction 

One of the most important physics issues at future hadron colliders, in particular at 
the CERN Large Hadron Collider LHC, will be the search for Higgs bosons. While 
in the standard electroweak theory (SM) one neutral Higgs boson is associated with 
electroweak symmetry breaking, the situation is more complex in many extensions 
of the SM. For instance, the well-known (non)supersymmetric two-Higgs doublet 
extensions of the SM predict the existence of three neutral spin-zero particles: one 
CP-odd and two CP-even Higgs bosons. The strategies to search for these particles 
[0] must take into account various possibilities. It may well be that the coupling to 
weak vector bosons of one of these Higgs particles is suppressed - or even zero, as 
in the case of a pseudoscalar - whereas its couplings to fermions is enhanced with 
respect to the SM Higgs couplings. 

In this article we investigate the signatures of a heavy Higgs boson in the ti 
decay mode at hadron colliders. Because of the large nonresonant ti background 
this is known to be a very difficult channel for the search of the SM Higgs boson. 
However, it is important to study this channel in a more general theoretical setting 
because, as mentioned above, Higgs particles with suppressed couplings to W and Z 
bosons but enhanced couplings to t quarks may exist. Such a scenario was analyzed 
in ||, and it was concluded in this paper that under favorable circumstances a 
scalar or pseudoscalar Higgs boson signal may be observable in the ti invariant- 
mass spectrum. Here we extend this analysis in several ways. We consider in the 
following a neutral Higgs boson (p of arbitrary CP nature - i.e., (p can be a scalar, 
pseudoscalar, or a state of undefined CP parity. We investigate in sect. 3 the signal 
of <p in the ti invariant-mass spectrum. We consider the case where this spectrum 
is reconstructed from the single lepton ti decay channel which is experimentally 
rather clean. Moreover, we exploit that in top quark production and decay the top 
spin is a "good observable" in the sense that due to the extremely short lifetime of 
the t quark the t spin-polarization and ti spin-spin correlations do not get severely 
diluted by hadronization and can be traced by appropriate correlations among the 
final state particles into which t and i decay [3-8]. The t spin-polarization and ti 
spin-spin correlations are sensitive to the Lorentz structure of the production vertex, 
and the decay ip — > ti leads to characteristic spin correlation signals ||. In sect. 4 
we analyze for the double-lepton ti decay channel the effects of (p production on ti 
spin-spin correlations, and we conclude in sect. 5. 

2 The ip —> ti channel including the QCD back- 
ground 

In order to be completely general we consider a neutral Higgs boson ip with unspec- 
ified CP parity. Its coupling to the top quark reads 
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where Gf is Fermi's constant, m t is the top mass, and a,a are the reduced scalar 
and pseudoscalar couplings, respectively. (For the SM Higgs boson a = 1, a = 0.) If 
act 7^ then ip has undefined CP parity, signaling CP violation in the Higgs sector. 



This can occur, for instance, in the two-Higgs doublet extension of the SM [JTO . 

In the following we consider a heavy p> with mass > 2m t with very small 
couplings to the weak vector bosons, but unsuppressed couplings to t quarks. Such 
a scenario may be realized in the supersymmetric extension of the SM || and also in 
nonsupersymmetric two-doublet models [l]. Then, for > 2ra* the mode <p —>■ ti 
determines the total <p decay width. To the extent that the decays of ip to W + W~ 
and ZZ are non-negligible the effects discussed in the next section become smaller. 

At high energy hadron colliders, for instance at the LHC, <p will be produced by 
gluon-gluon fusion through a virtual top quark loop |TT|. Scalar and pseudoscalar 
Higgs production by gluon fusion has been analyzed in detail in the literature ||12|| . 
The ip — > ti decay channel is affected by the large nonresonant tt background. 
The amplitudes gg — *> <p — > tt and gg — > tt interfere and produce at the parton 
level a characteristic peak-dip structure in the cross section [[IT], 0, [TjJ and in other 
observables if the tt invariant mass lies in the vicinity of the Higgs mass. As 
we are interested, apart from the tt invariant mass distribution, also in tt spin-spin 
correlation phenomena we have calculated the tt production density matrices for the 
reactions qq — > tt and gg — ► tt by taking into account the QCD Born diagrams and 
the s channel ip production diagram. The definitions and results are given in the 
appendix [T9I] . 



The tt spin-spin correlations are sensitive to the dynamics of top quark produc- 
tion. This may be exemplified by the the following examples. It is known [[22|, |3|, || 
that close to threshold the gluon-gluon fusion amplitude at Born level yields a ti 
pair in a ^So state. (The QCD Born amplitudes qq — > tt leave the ti pair in a 3 Si 
state.) Consider, on the other hand, gg — > <p — > ti. For a scalar and pseudoscalar 
p the ti pair is in a 3 Pq and in a 1 J S'o state, respectively. Let us now consider the 
expectation value of the product s t ■ sj of the t and i spins. A simple computa- 
tion yields that (s^-Sf) = —3/4 if ti is produced by gluon fusion at threshold, and 
(sfSf) = 1/4 (—3/4) if ti is produced by a (pseudo) scalar spin-zero boson and if the 
gg — > ti background is ignored. In Fig. 1 we have plotted, as a function of the ti 
invariant mass, this expectation value for gg — > ti without and with the inclusion 
of the s channel scalar, respectively pseudoscalar Higgs boson exchange diagram. 
The Higgs boson mass m v = 354 GeV was chosen to be close to the ti threshold 
of 350 GeV in order to exhibit the features discussed above. Above the threshold 
and resonance regions the Higgs effect in (s t -Sf) is actually more pronounced for a 
pseudoscalar than for a scalar, as can been seen from these figures. This will be 
further exemplified in sect. 4. 

The t and i quarks auto-analyze their spins by their parity- violating weak decays. 
We shall assume that t — ► W + b is the dominant decay mode, as predicted by the 
SM. It is well-known that the most efficient analyzer of the t spin is the charged 
lepton from subsequent W decay. Its spin analyzer quality is more than twice as 
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high as the W or b quark direction of flight. The corresponding t and i decay density 
matrices are given for instance in [23-25]. 

In the following section we consider two types of ti decay channels: first the 
"single lepton channel", where the t quark decays semileptonically and the i non- 
leptonically, 



t + t -> l + vib + qq'b 



(2) 



and vice versa. These channels have a good signature for top quark identification. 
Moreover, the t and i momenta can be reconstructed up to an ambiguity which can 
be resolved statistically ||26|| . These processes are therefore suitable for determining 
the ti invariant-mass spectrum. Second we consider the processes - which we shall 
call below "double lepton channel" - where both t and i decay semileptonically, 



t + t -> l + vib + l'~v v b . 



(3) 



As mentioned above, correlations among the directions of flight of the charged lep- 
tons are best suited to analyze ti spin correlations. We will therefore study angular 
correlations in the double lepton channel. 

For obvious reasons we demand the lepton in a semileptonic top decay to be 
either an electron or a muon. The SM predicts, to good approximation, a fraction 
of 24/81 of all ti events to decay into the single lepton channel and 4/81 into the 
double lepton channel. 



3 The tt invariant mass spectrum 

Throughout our calculation we use the narrow width approximation for ti produc- 
tion and decay, which is justified because of T t /m t , Y w /m w <C 1 and because we 
are primarily interested in normalized distributions and expectation values of ob- 
servables. Moreover we take the leptons and the light quarks including the b quark 
to be massless. Let us start with writing down our "master formula" for the cross 
section measure for the double lepton channel: 



da {pp -> ttX l + + ui + b + 1'~ + v v + b + X) 

M V [ 1 dx 1 [ dx 2 N x (x 1 )N- x (x 2 )e(s -Am 2 t ) 

xZT^gJO JO 
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In the first line on the r.h.s. Af = Bit — > 6/z/;) 2 is the square of the semileptonic 
t branching ratio, and N\(xi), Nx(x 2 ) are the parton distribution functions. The 
next four lines in eq. (Bj) represent the differential cross section at the parton level, 
written in terms of the t momentum direction k in the parton cm. system, the 
lepton momentum directions q ± , and the normalized lepton energies y± = 2E±/m t , 
defined in the t or i rest systems, respectively. The minimal value of the normalized 
lepton energies is \i — m^/m 2 . Forefactors are the strong coupling constant a s and 
the kinematic factor /3 = (1 — 4m;/s)^ 2 with the parton cm. energy s = XiX^s. 
The factor rj = (1 — yu) 2 (l + 2/x) /6 is chosen such that the lepton energy integrations 
are normalized to unity. 

The matrix element contains the coefficients A^ x \ B^ x \ . . . c$ X \ which depend 
only on s and the cosine of the angle between k and the beam direction, z = p-k. 
They are listed in the appendix. Note that only contributes to the rate, while 
the other coefficients lead to t and i spin polarization and spin-spin correlations. In 
particular, the total cross section for the parton subprocess AA — * tt is simply 

a^(s) = ^JdQ^ x \s,z), (5) 

which can be computed analytically. 

The last two lines in eq. (f|) finally give the distribution of b (b) energies E b (£5) 
and directions q b (q^), defined again in the t (t) rest system. The four 5-functions 
have their origin in the top and W propagators, for which we have used the narrow 
width approximation. One can see that once the lepton momenta are given, the b 
(b) momenta are fixed up to an azimuthal angle around the accompanying lepton, 
on which the matrix element does not depend. 

The cross section for pp collisions is simply obtained from eq. (|J) by substituting 
the parton distribution function iV^a^) with iVx(x2), which enhances the incoherent 
quark annihilation background. Note that eq. (|4]) can be as well applied to reactions 
where one or both of the t quarks decay hadronically; in this case the lepton and 
neutrino momenta are to be interpreted as the momenta of the (light) quarks into 
which the W decays, and one has to supplement the appropriate branching ratios. 



Let Mu = J(k t + kt) 2 be the invariant mass of the tt system. The spectrum 
da/dM t t is then obtained [@ from multiplying the total parton cross sections 
at s = M t | with the so-called luminosity functions 

L^\M th s) = ^ [ 1/T ^N x (rC)N- x ( T -) (6) 

where r = M t ij \fs- However, this approach does not account for the experimental 
difficulty of reconstructing the top momenta, which is present in the single lepton 
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channel, and it is not a priori clear to what extent a misidentification might disturb 
the signal. Apart from that, the implementation of cuts is not straightforward. 

We have therefore taken eq. (Q) as the starting point for writing a Monte Carlo 
generator. Importance sampling is performed in the variables x\ and X2, where we 
expect the biggest variations (due to the threshold behavior of the parton cross 
section and the form of the parton distribution functions). As can be seen from eq. 
(f|), the lepton energies can be simulated independently. Moreover, the construction 
of the b and b momenta requires only two (equally distributed) azimuthal angles, as 
mentioned above. 

The obstacle to reconstructing both the t and i momenta in the single lepton 
channel is the missing neutrino momentum. Four-momentum conservation and the 
various on-shell conditions determine the latter as the solution of a quadratic equa- 



tion up to a twofold ambiguity. In ref. a strategy was described how to improve 



on the odds for picking the correct solution: In a first step, cuts are applied on the 
rapidity of the top quark and the W boson on the hadronic side. Next one tests 
whether the reconstructed X\, X2 lie in the interval [0, 1]. Finally one chooses that 
neutrino momentum which - when combined with the lepton momentum - gives 
a value for the invariant l-v\ mass which is closest to the physical W mass. This 
procedure can raise the probability for a correct identification to 86-89%. 

It is a simple matter to implement this algorithm into our Monte Carlo generator: 
Although we know the correct neutrino momentum, we compute for each event the 
wrong solution and construct a "wrong" top momentum, which we select in 11-14% 
of all cases for a determination of M t j. In addition we apply the following cuts: 

Ht)\<3, p T (l)> 20 GeV, (7) 

for the rapidity of the top quark and the transverse momenta of the leptons. We 
will keep these cuts throughout below, also for the double lepton channel in the next 
section. For the numerics we put m t = 175 GeV and my/ = 80.4 GeV, and we set 
the partial widths T(ip — > W + W~ , ZZ) to zero. We took the parton distribution 
functions (PDF) from ref. [27], evaluated at the factorization scale A = m t , and 



convinced ourselves that our results do not change significantly if we vary A between 



mt/2 and 2m t or work with other PDF sets |£8 . 

Fig. 2 shows the resulting normalized invariant ti mass spectrum a~ l da / dM t i for 
pp collisions at a/s = 4 TeV in bins of 10 GeV, starting at the threshold of 2m t = 350 
GeV. In grey we show the result from the QCD background and superimpose as 
solid lines the curves which we obtained including the production of a Higgs boson 
of trip = 400 GeV. The respective values of the tpti couplings are chosen to be: 2a) 
a = 1, a = (scalar), 2b) a = 0, a = 1 (pseudoscalar) , and 2c) a = a = 1 (undefined 
CP parity). Each of the plots in Fig. 2 was produced with 50 000 events. This 
number results from assuming pp collisions with 10 fb _1 integrated luminosity and 
a detection efficiency of 1/3 for the single lepton channel. 

For the Yukawa couplings used in Fig. 2 the ip resonance becomes broader if 
a 7^ 0. As a consequence, the contribution to gg — > ti of a ip boson with pseudoscalar 
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component is, markedly above s > rrii, more visible than that of a pure scalar. The 
effect in the normalized spectra shown in Fig. 2 is therefore more pronounced for a 
if boson with a^O than for a pure scalar if. One can see that the typical peak-dip 
structure in the M t i spectrum described in |fjj survives at the level of the t and i 
decay products, even if one implements realistic cuts and misidentifies the top quark 
momentum in ~ 11% of all events. 

For the upgraded Tevatron an integrated luminosity of 10 fb _1 is eventually 
projected, but at ^/s = 2 TeV. At this cm. energy no significant signal of a if boson 
with the Yukawa couplings of Fig. 2 remains above the ti background. 

The situation clearly improves if we go to pp collisions at LHC energies, one 
reason being the suppression of the quark annihilation background. In Fig. 3 we 
show the same series of plots as in Fig. 2, but now for pp at 14 TeV. These plots 
were produced with 1.2xl0 6 events. This is a conservative number for the LHC with 
integrated luminosity of 100 fb~ . Qualitatively the signals are comparable to those 
in Fig. 2, but cleaner due to the larger number of events. The clear signals are of 
course due to the fact that we chose the Higgs mass m v = 400 GeV to be close to 
the threshold of 2m t = 350 GeV. To study the effect of different Higgs boson masses, 
we show in Fig. 4 a series of invariant mass distributions obtained for m v = 375, 
400, 450 and 500 GeV. In these plots, <f was taken to be pseudoscalar with Yukawa 
couplings a = and 5 = 1. 

Let us briefly discuss the statistical significance of our results, for the case of a 
pseudoscalar with mass 500 GeV (cf. Fig. 4d). Only 84.2% of the 1.2 x 10 6 events 
pass the cuts (|7|) and are filled into the histogram. If we consider the contents of 
the five bins below m^, i.e. the M t i interval between 450 GeV and 500 GeV, we 
find no =161 765 events for the background and 164 392 events if the if exchange 
is included, which is a 1.6% effect. Due to the large number of events this is a 
significant increase: With a 1 s.d. uncertainty of about ^/uq ~ 402, this translates 
into a 6.5 a effect. The statistical significance is, of course, higher for the curves 
depicted in Figs. 4a,b,c. 

The Higgs effect decreases if the decays if — > W + W~ , ZZ are non-negligible. If 
the couplings of if to the weak vector bosons are of SM strength then, using the 
same number of events, if masses, and Yukawa couplings as in Figs. 4a-d, the if 
resonance effect in the M t i spectrum becomes insignificant - except possibly for a if 
with pseudoscalar component and a mass close to 2m t . 

4 Spin-spin correlations of the ti system 

We have also studied the impact of the s-channel (^-exchange diagram on various 
correlations among the final state particles. We have concentrated on the double 
leptonic decay mode and investigated several observables involving the lepton ener- 
gies and momenta. In particular, the spin-spin correlation (s t -sj) discussed in sect. 2 
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suggests to investigate the observable 



= Q + Q_, (8) 

where Q ± are the unit vectors of the lepton momenta measured in the hadron cm. 
system. In other words, O is the cosine of the angle between the directions of flight 
of the two leptons as seen in the laboratory frame. Recall that the lepton momenta 
q ± in da ((4|) are defined in the t and i rest systems and thus related to Q± by two 
(rotation-free) boosts. 

The grey area in Fig. 5 shows the distribution of O for pp collisions at 6 TeV, 
as obtained from the background only. One can see that event topologies where 
the two leptons are back to back or parallel are slightly preferred. In the same 
figure we show as solid line the resulting distribution if we assume the presence of a 
pseudoscalar Higgs boson with mass m v = 400 GeV and Yukawa couplings a = 0, 
a = 2. In this case the distribution is shifted towards positive values of O. 

In order to quantify this effect one could study the asymmetry between the event 
numbers with positive and negative O. The same information is of course contained 
in the expectation value {O). The sensitivity of (O) to the couplings a and a is 



determined by the 1 s.d. width SO = y ((C 2 ) — (0) 2 )/N where N is the number of 
events in the sample. 

In Fig. 6 we show the theoretical expectation for (O), for pp collisions from 4 
to 14 TeV. The solid line and the dark grey band represent for a given cm. energy 
s/s the center value and the 3 a interval [(C) - 3SO, (O) + 350]. The width SO 
was computed with Nu = 2 x 10 5 events, which is 1/6 of the single lepton events 
used in the previous section. Also shown in Fig. 6 - as dashed line with light grey 
band - is the effect on (O) and SO caused by a = 400 GeV Higgs boson with 
pseudoscalar couplings a = 0, a = 2. One can see that this leads to an increase of 
(O) with a significance above the 3 a level. With the same coupling strength and 
mass, a scalar Higgs boson (i.e. a = 2, a = 0) does not yield a significant signal. 
This is expected from the discussion in sect. 2. 

Both Figs. 5 and 6 were produced with the cuts eq. (0). We observed that the 
signal gets enhanced by the cut on the lepton transverse momentum. 

One might think that the correlation effects increase if one works with the lepton 
directions q ± measured in the t and i rest systems. Although we consider the double 
lepton channel, it is possible to reconstruct both the t and i momenta 0, albeit with 
a certain loss in efficiency We have therefore investigated also the observable 

0' = q + -q_. (9) 

The results for (O) and (O') are qualitatively similar. However, in the case of (O') 
already the 2 a bands overlap, which means that the observable O' is less sensitive 
to the Yukawa couplings a and a than O. 

Finally we mention that neutral Higgs sector CP violation, i.e. ad ^ 0, can be 



tested with appropriate asymmetries and correlations p, |2lj, 25, 29| £3(J 
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5 Conclusions 



Extensions of the SM predict a number of neutral Higgs bosons tp. It is possible 
that one - or several - of these bosons have very small couplings to W and Z bosons 
but unsuppressed couplings to top quarks. In this paper we have considered the 
resonant production of such particles tp in the reaction p + p(p) — > t + t + X. We 
have studied the effect of tp on the tt invariant- mass distribution obtained through 
the single lepton channel. At the upgraded Tevatron (y/s = 2 TeV) with integrated 
luminosity of 10 fb _1 the signal of ip is insignificant - assuming that the Yukawa 
couplings of tp do not exceed SM strength - because of limited statistics and the 
large tt background from quark- ant iquark annihilation. For pp collisions at y/s = 
4 TeV, assuming again 10 fb~ x of integrated luminosity, tp bosons become visible 
in the tt invariant-mass spectrum for masses m v around 2 m t . For pp collisions at 
LHC energies and luminosities our analysis suggests that tp production is detectable 
for a range of masses m v . 

If such bosons should be found the analysis of their properties can be supple- 
mented by analyzing t and t spin-polarization and spin-spin correlation phenomena 
in the tp — > tt decay channel. The spin-spin correlation that we have studied is, at 
LHC energies, sensitive to a resonant tp contribution if tp has a sizeable pseudoscalar 
coupling to top quarks. This and other correlations and asymmetries should be 
useful in eventually pinning down the nature of these bosons. Our study was based 
on the lowest order parton matrix elements. Detailed analyses eventually require 
the inclusion of QCD radiative corrections, of hadronization and detector effects. 
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Appendix 



We define a decay density matrix p for the process t — > i + vi b in the narrow width 
approximation for the W boson by 



7T 



7Tl\yT\y 



S(Pw-m 2 w)p a > a = E (t*>\Tl\l + is £ b) (£ + u e b\T\t a ) , (A.l) 

His 



where pw, raw and IV are the four- momentum, mass and width of the intermediate 
W + , and a, a' are the t spin indices. The explicit form of p, conveniently evaluated 
in the top quark rest system, can be found e.g. in |25[ . The matrix p for the conjugate 
decay i — > l~i>i b is obtained from p by simple replacements. 

Similarly, we define the production density matrices for the partonic process 
AA — > ti (A = q, g) as 

4 A 2w = 4-^- E (t a i/3\T\XX) (\\\T^\t a >ip>) , (A.2) 

9s colors 

initial spins 

where the factor n^\) averages over spin and color of the initial state partons: 
= {2N c f = 36, n {g) = (2(iV£ - l)) 2 = 256. The factor g 4 s is taken out for 
convenience. A general discussion of the symmetry properties of these matrices and 
their decomposition in the t and i spin spaces is given in . Here we consider only 
QCD Born amplitudes and the lowest order gg — > ip — > ti amplitude. Therefore 
several of the coefficients in the general decomposition vanish, and the matrices 
are of the form 



(A.3) 



where a % are the Pauli matrices. In our case, the tensors C$ have the structure 
Qj = c o A) $ij + 4 A) £ ijl k + c[ X) pipj + 4 A) kikj + 4 A) {pikj + Pjk) ■ (A.4) 

Here k and p are the unit vectors of the momenta of the top quark and of the initial 
parton A, respectively, defined in the parton cm. system. (Neglecting transverse 
parton momenta implies that p is equal to the direction of the proton beam in 
the laboratory frame.) The squared matrix element of the reaction A A — > ti — ► 
l + uib + l'~viib is then proportional to 

Tr [pR {X) p\ = Pa'aR^l'wPP'P ' ( A - 5 ) 

which leads - after integration over the phase space and folding with parton distri- 
bution functions - to the representation of the differential cross section da given in 
eq. @. 
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The matrix gets only contributions from the QCD background, leading to 
the coefficients 





= N q 


2 
1- ^-(1 




2 v 


M 
c o 


= N q 


r P 2 

- -y' 1 


Si) 
c 4 


= N q 





cl 9) = N 



[3 2 + z 2 (l-xf 



c { 6 q) = N q [-z(l-x)} 
= 4 9) = 0, 



(A.6) 



with 



x 



2m t 



and the normalization factor 



P = VT 



X" 



z = P k , 



TV = — — ^ 
9 4 • 



(A.7) 
(A.8) 



As usual, the SU(N C ) color factors are C F = (N 2 , - 1)/(2N C ), C A = N c . 

We split the coefficients of the gluon fusion density matrix into contributions 
from the QCD background and contributions involving the if exchange diagram: 



A(g) — a(9) I Ma) 

71 _ ^Born + ^-w ' 



Born 1 " "ip 

and similar for the other coefficients. For the background part we obtain 



(A.9) 



i Born 



c 0,Born — 



J9) 

"4, Born 



-"5, Born 



JV(0) 



2^2 



2 1-/3^ 
1 + x 2 ^~ 



X 



1 - (3 2 z 2 



X 



1 - f3 2 z 2 



X 



1 - (3 2 z 2 



1 - /5 2 ^ 2 
1 — X 



9 1 + x 



_ n , „ 2 2 x 2 (3 + 4x + 2:r 2 ) 

1 + 2x(l + x) + /3V - 



1 - /3V 



(A.10) 



c 6,Born 



B 



(<?) 
Born 



= N®[-z(l-x)] 



c 2,Born — u ■ 



P 2 (l-z 2 ) 



1 - (3 2 z 2 



For the contributions from if exchange it is convenient to define the following two 
functions, 



D = (-2 + (3 2 sC 



D = sC 



s - m 2 + im v T v ' 



(A.ll) 



s -m 2 + im^Tp ' 
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which contain the loop integral 



Co = Co(p 1 ,p 2 ,mt,m t ,m t ) = — 

2s 



In 



.1-/3, 



(A.12) 



as well as the Breit-Wigner form of the ip propagator. With the help of D and D 
the remaining contributions can be written in a compact way: 



la) 

c 0,<fi 



M 
la) 

c 5,<p 

la) 

c 4,¥> 



= N^ 2 
a 

+ 

= N? 

= Nf 

+ 
= 

= 4" 

= la) 



(3 2 a 2 Re{D) + a 2 Re{D) 



(0 2 a 2 + a 2 ) [a 2 \D\ 2 + ~a 2 \D\ 
((3aa)lm(D - D) , 
~(3 2 a 2 Re(D) - a 2 Re(D) 
{(3 2 a 2 -~a 2 ) (a 2 |D| 2 + a 2 |D| 
(-/3aa)Re(D + D) + Nf\-2l3aa) (a 2 \D\ 2 + a 2 \D\ 
(-2/?V)Re(L>) + N^\-2(3 2 a 2 ) 



a z \D\' z + a 2 \D 



. 



(A. 13) 



Here a and a are the reduced scalar and pseudoscalar Yukawa couplings defined in 
eq. ([!]). The normalization factors are 



C A l-f3 2 z 2 4C F/ 

m t \ 2 

v J 1 

i 



2 C F C A 32tt 2 
11/1 



X 



(5 2 z 2 



(A.14) 



v 



X 



2C F V32tt 2 , 

with v = {^Gp)- 1 / 2 = 246 GeV. The coefficients given in this appendix agree with 
those obtained in B. The cross sections for the partonic reactions are determined 
by A^ x \ and our expression for da(gg — ► tt)/dz agrees with [[|. 
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Figure 1 




Fig. 1: Expectation value (s t • sj) of the product of t and i spins in the process 
gg — > ti, plotted as function of the parton cm. energy V§. The solid line 
is the result for the background, while the dashed (dotted) line shows the 
effect of a scalar (pseudoscalar) Higgs boson with mass = 354 GeV 
and reduced Yukawa couplings a — 1, a — (a — 0, a — 1). 
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Figure 2 
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Fig. 2: Normalized tt invariant mass 
spectrum for pp collisions at 4 TeV. 
Shown is the background contribution 
(grey) and the effects of a Higgs boson 
with mass m v = 400 GeV: a) a—1, a = 
(scalar), b) a = 0, a=l (pseudoscalar), 
c) a — a — 1 (mixed). 
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Figure 3 




Figure 4 




Fig. 4: Normalized tt invariant mass spectrum for pp collisions at 14 TeV. Shown 
is the background contribution (grey) and the effects of a pseudoscalar 
Higgs boson with couplings a — 0, a — 1. The different plots a), b), c), 
d) are for mass values m v = 375, 400, 450, 500 GeV, respectively. 
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Figure 5 




Fig. 5: Normalized distribution a -1 da / dO , where O is the cosine of the angle 
between the two leptons in the double lepton channel of pp collisions at 
yHi = 6 TeV. The grey area is the contribution from the background only. 
The solid line shows the effect of a pseudoscalar Higgs boson with mass 
m v = 400 GeV and couplings a = 0, a = 2. 
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Figure 6 
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Fig. 6: Expectation value (O) of the cosine of the angle spanned by the two leptons 
in the double lepton channel of pp collisions, plotted as function of the cm. 
energy y/s. The solid line and the dark grey band represent the expectation 
from the background only, for 2 x 10 5 events at the 3 a level. The dashed 
line with the light grey band was obtained with a pseudoscalar (a = 0, 
a — 2) Higgs boson with mass = 400 GeV. 



19 



